A model-independent analysis of the dependence of 
the anomalous J/ip suppression on the number of 
participant nucleons 



A. P. Kostyuk 1 2 and W. Greiner 1 

1 Institut fur Theoretische Physik, Universitat Frankfurt, Germany 

2 Bogolyubov Institute for Theoretical Physics, Kyiv, Ukraine 

Abstract. A recently published experimental dependence of the J/ip to Drell-Yan ratio on 
the measured, by a zero degree calorimeter, forward energy Ezdc m Pb+Pb collisions at 
the CERN SPS is analyzed. Using a model-independent approach, it is shown that the data 
are at variance with a earlier published experimental dependence of the same quantity on the 
transverse energy of neutral hadrons The discrepancy is related to a moderate centrality 
region: 100 < N p < 200 (N p is the number of participant nucleons) and is peculiar only to the 
data obtained within the 'minimum bias' analysis (using the 'theoretical Drell-Yan'). This 
could result from systematic experimental errors in the minimum bias sample. A possible 
source of the errors is discussed. 



Recently, the NA50 collaboration published new data on the centrality 
dependence of the J /ip suppression pattern in Pb+Pb collisions at CERN SPS 1 . 
The centrality of the collisions was estimated by measuring the energy of projectile 
spectator nucleons Ezdc by a zero degree calorimeter. It was, however, mentioned 
|2] that the new data may be at variance with the ones published by the same 
collaboration in Refs. where the transverse energy of produced neutral 

hadrons Et was used as a centrality estimator. The two sets of data (for brevity, 
we shall call them respectively 'Ezdc data' and 'Et data'), having very similar 
qualitative behavior, seem, nevertheless, to be at variance quantitatively. It was 
demonstrated in Ref. [2] that using three completely different J/ip production models 
if a model is fitted to the Et data, it does not agree with the i?z£>c-that and vice 
versa. Earlier it was mentioned joj that the co-mover model, which describes very well 
the Et data, does not show similar agreement with the Ezdc se t- 

Although the previous analysis pointed out a possible problem, its results can still 
be considered as inconclusive. Indeed, one may believe that the problem is related to 
the used models rather than to the data, and the 'truly correct model' should agree 
with the both data sets simultaneously. 

In this letter, we use a model-independent analysis of the data to show that any 
model would be at variance at least with one of the data sets. This points out presence 
of systematic errors. 

There are two types of data in each set, depending on the analysis procedure used 
to obtain the ratio R of J/ip multiplicity to the multiplicity of Drell-Yan pairs. In the 
standard analysis, the number of both J/ijj and Drell-Yan events are extracted from 
the measured dimuon spectra, i.e. the physical Drell-Yan is used as a reference. In 
the 'minimum bias' analysis, only the number of J/ip events and the total number of 
collisions (minimum bias events) in each centrality bin are found from the experiment. 
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The ratio R is obtained by dividing the measured ratio of the number of J j ip events to 
the number of minimum bias events by the theoretical ratio of Drell-Yan multiplicity 
to the probability of a minimum bias collision at the corresponding centrality. Then 
the overall factor is fixed using the standard analysis data. The theoretical Drell-Yan 
to minimum bias ratio is given by 

inY] _ I d 2 bP(E x \b) {DY(b))P int (b) 

[ > E * J <PbP{E x \b)P int {b) ' U > 

where (DY(b)) is the average number of Drell-Yan dimuon pairs per Pb+Pb collision 
at impact parameter b and Pint (b) is the probability that an interaction between two 
nuclei at impact parameter b takes place. Both quantities are calculated in the Glauber 
approach. The subscript x stands for T or ZDC. 

The ^^-distribution of events at fixed impact parameter b, P(Ex\b), is assumed 
to have a Gaussian form with the central value and dispersion given by 



(E T (b)) = qN p (b) a ET (b) = q^aN p (b). (2) 

Here N p (b) is the average number of participant nucleons in a Pb+Pb collision at 
impact parameter b calculated in the Glauber approach. The parameter values 
q = 0.274 GeV and a — 1.27 [7] are fixed from the minimum bias transverse energy 
distribution [5]. 

The i?zDC-distribution P{EzDc\b) is assumed to be a Gaussian with the central 
value and dispersion given by Eqs. (1) and (2) of Ref. pQ. 

At a given dependence (J/ip(b)) of the average number of J ftp's per Pb+Pb 
collision on the impact parameter b, the ratio R in a finite-size E x -bin, £™ m < E x < 
E™ ax can be calculated as: 

/5T dE x J d 2 bP(E x \b)(J/iP(b))P mt (b) 

R(E™ m <E X < E™ ax ) = i ( 3 ) 

Jj in dE x jd?bP{E x \b){DY{b))P int {b) 

In absence of essential systematic errors, the Et- and Ezdc data should be consistent 
with each other: the formula © should fit the both data sets simultaneously: the % 2 
per degree of freedom, x 2 fndf, that includes Et as well as Ezdc data points should 
be in the vicinity of 1. In the following we estimate in a model- independent way the 
lower bound for x'/ndf and show that it is in fact essentially larger than 1. 

Instead of x 2 P er degree of freedom, we shall calculate the x 2 per experimental 
point, x 2 fnep, which provides a lower estimate for \ 2 f n df- 

X 2 /ndf > x 2 /nep, (4) 

because ndf < nep. 

We approximate the unknown function {J/ip(b)} in the following way: 

• We fix J nodes bj on the &-axes: bo = 0, while bj is taken sufficiently large 
(6,7 = 22 fm) so that the probability of interaction of two Pb-nuclei and nuclear 
modification of J/ip production at b = bj are negligibly small. In our calculations, 
we fix the rest of nodes in two different ways: 

(a) equidistant 

bj+i — bj = Ab = const, j = 0, . . . , J — 1 (5) 

and 

(b) so that the difference of the average number of participant nucleons between 
every two neighboring nodes is the same: 

N p (bj) - N p (b j+ i) = AN P = const, j = 0, . . . , J - 1. (6) 
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• The values of (J/ip(b)) in all nodes except the last one: 

(J/WM^xj, i = 0,...,J-l, (7) 

are treated as free parameters. The value in the node bj is fixed assuming the 
same J /if) to Drell-Yan ratio in extremely peripheral Pb+Pb collisions as in p+p: 

(JMM> = *J = 53-5 (DY(bj)) (8) 

• The value of (J/ip(b)) between the nodes is given by a linear interpolation: 

(j/m) = {bi+1 - b l Xi + {b h - bi)Xi+1 , b j<b <b j+1 . ( 9 ) 

bj+i- bj 

• There is no reason to expect oscillating dependence of the J/ip multiplicity on 
the centrality. Therefore, Xj QiTG subjected to the constraint 

< ?Z (10) 

{DY{b 3 ^)) ~ {DYib^Y 1 > 

which ensures monotonicity of the J/ -0 to Drell-Yan ratio R as function of b. 

• One more constraint is needed to avoid negative values of J/ip multiplicity: 

x >0. (11) 

Any non-negative function (J/ip(b)) related to a monotonic R(b), including even a 
disconnected one, can be arbitrarily well approximated in the above way, provided that 
J is sufficiently large. Therefore, minimizing of x 2 with respect to all Xj, < j < J, at 
sufficiently large J one gets the value, which is, at least, not larger than the minimum 
X 2 in any reasonable model. 

Using Eqs.(JHJ and one gets a linear dependence of R on xy. 

R (n) =Y,A^ Xj . (12) 

3=0 



Here 



, , JpM.min dE x a^(E x ) 

4 n) = — — ^ ; (13) 

dE x / °° dbbP(E x \b){DY(b))P mt (b) 



S^SZZ dE x [a-{E x ) + a+(E x )} 

dE x r °° dbbP(E x \b)(DY(b))P mt (b) 



4 n) = — : — — . i = i,-. .,•/-!; (14) 



where 



jftSrr dE x aj(E x ) 

4 n) - p(B , mM 54 ; (is) 

J5w.».» dS« / °° dbbP{E x \b){DY{b))P int (b) 

a+(E x )= [ 3+1 dbbP{E x \b) hl+1 ~, 6 F mt (6), j = 0,...,J-l; (16) 
a7(£*) = dbbP(E x \b)^—^f--P mt (b), j = l,..., J. (17) 

A,-! - ! 

In our fitting procedure, the initial values of Xj 9XC chosen randomly (using a 
random number generator) . Then the x 2 is minimized in such a way that the constraint 
ifTUjl is fulfilled at each iteration step. 
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The result of the fit of the NA50 minimum bias data is presented in Tab. D and in 
Figs^andl^l As is seen from the table at J > 100 the value of x 2 / ne P saturates 
and ceases to depend on the way, how the nodes are fixed. Therefore J = 200 
can be accepted as sufficiently large to provide a reliable lower estimate for \ 2 /nep. 
First, we fitted the Et- and Ezdc data separately ( 1 Et fit' and the 1 Ezdc fit', 
respectively). In both cases the lower estimate of x 2 /nep is smaller than 1, which 
indicates that the difference between the data from 1996 and 1998 runs within each 
minimum bias dataset can be well explained by statistical fluctuations only. But if one 
compares the Et data and Et dependence of the J/ip to Drell-Yan ratio, calculated 
with Xj's found from the Ezdc fit (see Fig. 0), one can observe clear discrepancy 
at moderate centrality (N p — 100-200). Similar discrepancy is observed, when the 
Ezdc dependence resulting from the Et fit is compared with the Ezdc data. 

Finally, we tried to fit the both sets simultaneously (the 'combined fit'). The 
calculations show that x 2 /nep = 3.28. This means that y/ 2 /ndf > 3.28 in any model. 
Again, it is clearly seen from Figsn an d[21that the problem is related to the centrality 
region N p = 100-200. The discrepancy cannot be explained by statistical fluctuations 
(this hypothesis is excluded at the level of about 10 -19 !). Therefore, at least one of 
the data sets contains essential systematic errors. 

There are reasons to suspect that these are the Ezdc data that contain the 
systematic errors and the source of these errors is a distortion of the minimum bias 
sample. Indeed, a crucial point of the minimum bias analysis is the assumption that 
both the measured Drell-Yan multiplicity as well as the experimental minimum bias 
centrality distribution are exactly proportional to the corresponding theoretical values, 
i.e. to the numerator and denominator of Eq.([Q), respectively. The experimental and 
theoretical minimum bias probabilities should cancel (up to a constant factor), when 
one divides the experimental J /i\> to minimum bias ratio by the theoretical Drell-Yan 
to minimum bias ratio leaving finally only J/ip to Drell-Yan ratio. 

While the experimental minimum bias Et distribution perfectly agrees with the 
Glauber model (see Fig. 1 of Ref. jSI ) :t the agreement of the Ezdc minimum bias 
distribution (see Fig.l of Ref. U) seems to be not so perfect. Unfortunately, the 
data plotted in Ref. 1 have never been published in numerical form and the value of 
X 2 /dof characterizing the quality of the fit of the minimum bias Ezdc distribution 
has not been quoted. Nevertheless, it is seen from the plot that the experimental 
points lie slightly above the theoretical curves in the region 25 > Ezdc ^ 18 TeV 
(which corresponds to 100 < N p < 200). Although the difference does not exceed even 
the size of the point symbols, it, nevertheless, indicates a sizable discrepancy, due to 
logarithmic scale of the vertical axis of the plot. Because of this discrepancy, the 
measured and calculated minimum bias distribution do not cancel completely, which 
may lead to the observed systematic errors. 

The question about the origin of the distortion of the minimum bias distribution 
is not simple. The final answer can hopefully be given after a detailed analysis of the 
raw experimental data supplemented by Monte-Carlo simulations. This is out of the 
scope of our present analysis. We restrict ourselves only to formulating a 'working 
hypothesis', which is to be tested in course of more detailed investigations. The 
observed discrepancy may come from interactions of reaction fragments with nuclei 
inside the zero degree calorimeter. If such an interaction takes place, a fraction of the 

| Disagreement is observed only in the low Et region, because the efficiency of the target 
identification algorithm is less than unity there. 
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produced particles (those having large transverse momenta) may leave the calorimeter 
without deposing their energy in it. This may cause losses of the zero-degree energy, 
which can result in distortion of the shape of the minimum bias distribution. If our 
hypothesis is correct, the data analysis procedure should be modified to take into 
account the losses of the zero degree energy. Otherwise the mentioned problem will 
persist even in a minimum bias analysis of the new, year 2000 data. 

From Fig. one can easily see that the minimum bias sample is much more 
sensitive to this effect than the dimuon one. Indeed, if spectators from a peripheral 
collision loose an essential part of their energy, this event is registered by the apparatus 
as a moderate centrality collision. The number of peripheral minimum bias events is 
about an order of magnitude larger than the number of events at moderate centrality. 
Therefore, even if the fraction of peripheral events that suffer from losses of the zero 
degree energy is small, the contamination at moderate centrality can be sizable. In 
contrast, the number of Drell-Yan events at moderate centrality is larger than the 
number of peripheral ones. Therefore the Drell-Yan sample is much less influenced. 
This is supported by our comparison of the standard analysis data. Even the newest 
data which have smaller statistical errors than older ones, do not show any 
contradiction between Et- and the Ezdc data sets: according to our analysis, the 
lower bound of % 2 /ndf for the combined fit of the new standard analysis data is 0.62. 

The region of low Ezdc {Ezdc 9 TeV), i.e. the region of the 'second drop' 
of the Ezdc data, deserves special attention. All models that have been successful 
describing the E T data at E T > 100 GeV 01113111111121 explain the 'second drop' of the 
J/ip to Drell-Yan ratio by .E^-losses in the dimuon event sample with respect to the 
minimum bias one and/or by the influence of -&r-fluctuations on the J/tp production 
or suppression. Both effects do not influence the Ezdc data, this means that all 
the models 0110111111121 fail to describe the E Z dc data at E Z dc < 9 TeV. If the 
drop at low Ezdc indeed exists, the correct model should explain the high-centrality 
behavior of both Et- and Ezdc data by a sharp decrease of the J/ip multiplicity at 
small impact parameter. 

It is, however, not unlikely that the 'second drop' in the Ezdc data is an artifact of 
the data distortion. Indeed, Fig. 1 of Ref . 1 shows clear disagreement of the minimum 
bias distribution with the Glauber model at low Ezdc- This, as was explained 
above, invalidates the minimum bias analysis procedure. Therefore, the apparent 
agreement between the Et- and Ezdc data in high centrality region obtained within 
our analysis (see Figs and |2J may be purely accidental. In this case, there would be 
no contradiction with the models El HH Q2| ■ Hence, refining of the Ez dc data is 
very important. It would be hopefully able to narrow the set of successful models in 
this field, which recently became too broad. 

In conclusion, we performed a model independent comparison of two data sets 
published by the NA50 collaboration: the Et- and E'zDC-dependence 01 the J/ip to 
Drell-Yan ratio R in Pb+Pb collisions at SPS. We have shown that the two sets are 
at variance. This problem seems to be related only to the data obtained within the 
minimum bias procedure. A possible reason for the problem may be a distortion of 
the minimum bias data sample. We argue that refining of the data is important to 
clarify the qualitative and quantitative behavior of the ratio R at Ezdc ^ 9 TeV. 
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Figure 1. The dependence of the J/ip to Drell-Yan ratio R on the transverse energy Et- The 
corresponding average number of participant nucleons N p are shown on the upper axis. The 
points with error bars are the NA50 'minimum bias' analysis data (the Et data). The dotted 
and dot-dashed line are fitted to the Et data and the Ezdc data (see Fig. |5J, respectively. 
The solid line represents the best possible fit to the both data sets simultaneously. 
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Figure 2. The dependence of the J/i/> to Drell-Yan ratio R on the energy of zero degree 
calorimeter Ezdc- The corresponding average number of participant nucleons N p are shown 
on the upper axis. The points with error bars are the NA50 'minimum bias' analysis data 
(the Ezdc data). The meaning of lines is the same as in Fig. 
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Figure 3. The probability distribution of minimum bias and Drell-Yan events as a function 
of the number of spectator nucleons N sp in arbitrary units (a.u.), obtained within the Glauber 
approach. 



A model-independent analysis. 



10 



Table 1. The value of \ 2 P er experimental point at different values of the number of nodes 
J and for two different ways of choosing the nodes (a) and (b) (see text for details). 



J 




fit 


Ezdc fit 


Combined fit 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


3 


6.62 


3.93 


21.6 


69.9 


25.5 


51.8 


6 


2.07 


1.60 


1.31 


6.19 


3.96 


6.72 


12 


0.84 


0.83 


0.86 


1.14 


3.44 


3.46 


25 


0.69 


0.68 


0.82 


0.83 


3.32 


3.34 


50 


0.67 


0.67 


0.81 


0.81 


3.29 


3.29 


100 


0.67 


0.67 


0.80 


0.80 


3.28 


3.29 


200 


0.67 


0.66 


0.80 


0.80 


3.28 


3.28 



